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We theoretically and experimentally generate stationary crescent surface solitons pinged to the 
boundary of a micro-structured vertical cavity surface emission laser by using the intrinsic cavity 
mode as a background potential. Instead of a direct transition from linear to nonlinear cavity modes, 
we demonstrate the existence of a symmetry-breaking crescent waves without any analogs in the 
linear limit. Our results provide an alternative and general method to control lasing characteristics 
as well as to study optical surface waves. 
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By state-of-the-art semiconductor technologies, micro- 
cavities have provided a controllable confinement and 
manipulation of photons with small mode volumes and 
ultrahigh quality factors 0, • For an integrable cav- 
ity shape, supported stable resonance modes have at- 
tracted much attention in photonics, quantum electrody- 
namics, and telecommunications, due to their potential 
applications to modulate spontaneous emission and make 
thresholdless lasing [3|. If the cavity shape is slightly 
deformed into a chaotic one, unstable localized waves, 
coined as scar modes, are found also to support lasing, 
which produce unidirectional outputs and provide an al- 
ternative understanding in the correspondence between 
classical and quantum system 0-0] ■ 

To form standing light waves, optical cavities are typi- 
cally defined by the geometry of media with a higher re- 
fractive index. Confined by total internal reflection at the 
interface, whispering-gallery modes arc almost grazing 
incidence patterns, which act as filters, delay lines, cou- 
plers, and sensors for a broad areas from optical commu- 
nications, information processing, to biophotonics [H, 
Instead of the cavity modes with periodic orbits, light 
wave can attach to the boundary of materials under ap- 
propriate conditions [Iol - [T2j ]. As electrons localized at 
crystalline surfaces known as Tamm and Shockley states 
[HI IHl' optical surface waves are localized at the in- 
terface between two different media. Direct observa- 
tions of optical surface states have been demonstrated 
in photonic lattice edges [ll| and periodic waveguide ar- 
rays [16j . 

Without introducing any symmetry-breaking in the 
geometry, optical surface waves can also be induced 
through optical nonlinearity, resulting in the formation 
of surface solitons [13 • These surface solitons are sup- 
ported even in uniform materials, and can not find any 
analogs in the linear cases. A specific type of localized 
states, in the form of crescent surface solitons, is pinged 
to a circular boundary with a similar shape as Barchan 
sand dunes [l8[ . Optical crescent solitons have been 
dieted theoretically to exist in highly nonlocal media 
through the superposition of two vortex beams with dif- 
ferent topological charges (2(j, or by introducing inho- 
mogeneous losses [2l|. In rotating Bose-Einstein con- 
densates, matter- wave crescent vortex solitons are found 
in the 2D Gross-Pitaevskii equation combining the lo- 
cal self-attractive nonlinearity and a quadratic-quartic 



potential [22|. With a structure of concentric rings in 
the refractive index modulation, stationary and rotating 
surface solitons with thresholdless formation powers are 
shown to exist at the edge of guiding structures [23j . 

In this Letter we demonstrate a new type of micro- 
cavities, fabricated on the surface of a Vertical Cavity 
Surface Emission Laser (VCSEL), and report experimen- 
tally the formation of crescent surface solitons by col- 
lecting near-field radiation intensities. By increasing the 
injection current, we analyze transitions between linear 
Laguerre-Gaussian-like cavity modes and nonlinear opti- 
cal patterns resembling soliton rings . In contrast to a 
direct crossover between linear and nonlinear modes, we 
introduce the concept to design the intrinsic linear cavity 
modes as a background potential, which near the thresh- 
old lasing condition supports stationary single-, double- 
, and quadruple-humped crescent solitons without any 
counterparts in the linear limit. Our numerical results 
based on a nonlinear wave equation within a focusing 
medium are in a good agreement with the experimental 
observations. The experimental and numerical investiga- 
tions in this work provide an effective state for investi- 
gating surface waves in micro-structured semiconductor 
lasers. 




FIG. 1: (Color online) (a) Schematic diagram of the device 
structure, (b) The L-I curve, power versus current, for our 
VCSEL with a surface microstructure. The inset shows the 
top- view SEM image, with the geometric diameter D = 2fim, 
two radius R2 = 6.5/xm and R3 — 7.5/im. 

Our optical cavity is fabricated on an ion-implanted 
VCSEL, with the schematic diagram shown in Fig. 1(a). 
The epitaxial layers of the VCSELs are grown by metal 
organic chemical vapor deposition (MOCVD) on a n + - 
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FIG. 3: (Color online) (a) Experimental demonstration of 
the mode pattern resembling a soliton-ring structure at the 
injection current 35mA. (b) shows the numerical result for a 
4 th -order nonlinear cylindrically symmetric mode composited 
by the superposition of unstable eigen-modes with azimuthal 
indexes of (d) I = 3, (e) I = 7, (f)Z = 11, and (g) I = 14, re- 
spectively, (c) shows the corresponding azimuthal instability 
spectrum. 



FIG. 2: (Color online) Experimental demonstration of a mode 
transition for crescent surface waves with the injection cur- 
rents for (a) 3mA, (b) 5mA, (d) 6mA, (e) 7mA, and (f) 
9.5mA, respectively, (c) shows the corresponding numerical 
result for the supported linear LGij-like cavity mode. Non- 
linear crescent surface modes with (g) single-, (h) double-, 
and (i) quadruple-humps are shown with the same simulation 
parameters but just varying the injection current. 



GaAs substrate, with a graded-index separate confine- 
ment heterostructure (GR1NSCH) active region formed 
by undoped triple-GaAs-AlGaAs quantum wells placed 
in one lambda cavity, where the emitting window is de- 
signed as 15/im in diameter and the emitting aperture is 
confined by the implanted region as 20//m in diameter. 
On the surface, we etch the emitting window by the fo- 
cus ion beam (FIB) to lfim in depth with a donut-shape 
mask of lfim inner radius and 6.5^m outer radius. From 
the scanning electron microscope (SEM) image, the inset 
in Fig. 1(b), top view of the surface structure displays 
that in the center there is a circular region with a higher 
refractive index, surrounded by the etched donut-shape 
medium with a lower refractive index, and then comes the 
outer annular cavity with a higher refractive index mod- 
ulation again. The L-I curve, light versus current, of the 
surface-structured VCSEL is shown in Fig. 1(b). The 
threshold current of the device is about 6mA with the 
lasing wavelength at 855nm. It would be demonstrated 
later, this kind of cavity geometry can generate crescent 
surface waves pinging to the outer circular boundary ef- 
fectively 

Next, we measure the near-field electromagnetic inten- 
sity distribution at a fixed injection current by a charge- 
coupled device (CCD) camera through a standard micro- 
scope with a 100A lens. While the VCSEL is operated 
below threshold, for example, at the current of 3mA, it 



can be seen clearly in Fig. 2(a) that spontaneous emission 
pattern just reflects the geometry of our optical cavity in 
the VCSEL. Since the circular cavity in the central region 
has a smaller mode volume and a higher quality factor, 
as the injection current increases, it is clearly to see that 
the light intensity in the central region becomes brighter 
and brighter, as shown in Fig. 2(b). After the VCSEL is 
turned on, our device begins to lase. Instead of a single 
central spot, a series of surface wave patterns resembling 
crescent surface modes are presented in Figs. Hfd-f) for 
the injection currents of 6, 7 and 9.5mA, respectively. 
Even though our optical cavity is designed with a sym- 
metric circular geometry, these localized surface modes 
have a symmetry-breaking shape along the azimuthal di- 
rection. Moreover, these lasing modes are generated at 
the outer circular boundary first, then multiple-humped 
profiles in the shape of crescent waves appear regularly 
in the etched donut-shape region. When the operation 
current is selected well above the threshold, nonlinearity 
becomes a dominant effect in the formation of the trans- 
verse optical patterns, such as 35mA shown in Fig. GOa). 
The observed structure itself is believed to resemble a 
cluster of soliton rings 24]. This is confirmed by the 
subsequent theoretical analysis summarized below. 

The complex spatiotemporal dynamics and the pat- 
tern formation in a broad-area semiconductor laser cav- 
ity have been theoretically analyzed to study filamentary 
behaviors, soliton manipulations, and transverse mode 
transitions. To illustrate theoretically the formation of 
observed crescent surface modes in a semiconductor laser 
cavity, we start with the coupled equations for the slowly 
varying complex envelope of the electric field E and car- 
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rier density N [251 ] : 

d t E = -{l+n + iO)E - i2CB(N -l)E + i d x 2 E, 
d t N = -7 [N + /3N 2 - I + \E\ 2 (N - 1) - dd ± 2 N] , 

where C is the saturable absorption coefficient scaled to 
the resonator transmission, d± 2 is the transverse Lapla- 
cian describing the diffraction in the paraxial approxi- 
mation, rj is the linear absorption coefficient due to the 
material in the regions between the semiconductor and 
the mirrors, 9 is the cavity detuning, — (i + a) rep- 
resents the absorptive response of material, a is the cor- 
responding line-width enhancement factor, and d is the 
diffusion constant of the carrier scaled to the diffraction 
coefficient. 7 and (3 are the normalized decay rates of 
the carrier density that describe the non-radiative and 
radiative carrier recombinations, respectively. External 
injection current is denoted by I. 

To find spatially localized solutions numerically, we as- 
sume that the diffusion length of carriers is much smaller 
than the diffraction length of the electromagnetic field, 
i.e., d <§C 1, also neglecting the second-order effects in the 
evolution of the carrier density by setting /3 = 0. When 
the system reaches equilibrium, as it happens in our laser 
cavity, the field \E\ 2 becomes stationary and we can re- 
place N — 1 to the first order in Fourier space [26j], i.e., 



N -1 



cxp [-7(1 + \E\ 2 ~ dd 2 )(t - t')] 7 (I - l)dt'. 



By neglecting the variation in the coefficient of differen- 
tial operator and assuming d constant, we derive a re- 
duced dissipative wave equation for the electromagnetic 
waves in the semiconductor microcavity p7j . 



id T E + 66 E — (1 + iS) 



-(1 + 77) + 



2C(J- 1) 
1 + \E\ 2 



+ (1 - iSd)d 2 ± E + V{r, 6)E = 0, 



E 



(1) 



where r = at and the spatial coordinates is normal- 
ized with a factor We also define 9 = 9 — 9q for 
#0 + ce(l + rj) = 0, and 6 = 1/a to simplify the notations. 
The last term in Eq. JT]), V(r, 0), is added to account 
the refractive index modulation for the cavity geometry. 
When the line- width enhancement factor a ^> 1 (6 <C 1), 
the proposed equation (jTJ) can be approximated by the 
generalized nonlinear Schrodinger equation, where mod- 
ulation instability is known to lead to the formation of 
nonlinear patterns and solitons (28j . 

Before applying the nonlinear wave model, we define 
the index modulation based on the lateral geometry mea- 
sured experimental in Fig. Ufa) by assuming the effective 
refractive indexes as 3.49 and 2 for the regions without 
and with the FIB etch process, respectively. The cal- 
culated linear eigenmodes (TE modes) of this cavity ge- 
ometry with the lasing wavelength are found to support 
two kinds of cavity modes; one has a 4-rings profile very 
similar to a Laguerre-Gaussian(LG) mode of 4*' l -order in 
the radial and th in the azimuthal directions, denoted 



by LG 1 ^ , for which the field distribution occupies most 
area of the whole cavity, as shown in Fig. (He); while 
the other one has a single peak just inside the central 
region that is close to a LGq. These modes are typical 
eigenmodes in the optical fiber waveguides and circular 
resonators. 

The central region in this micro-structured geometry 
acts as an auxiliary cavity, of which the size plays a cru- 
cial role in our design to generate crescent surface modes. 
In the limit of a geometry without this central cavity, our 
VCSEL is a broad-area circular cavity, with an annual in- 
dex modulation in the outer ring. Although this type of a 
broad-area cavity also supports radially distributed lin- 
ear modes as LG^-like mode, but practically it is hard 
to produce a symmetry-breaking mode as crescent waves 
due to the mode competitions. On the other hand, if 
the size of the central region is larger, we found experi- 
mentally, not shown here, that only lasing in the central 
cavity is possible and the lasing in the outer region is 
totally suppressed. 

Due to the fact that the intrinsic LG^-like mode in 
Fig- is suppressed as a result of mode competition 
with the single-peak mode in Fig.^a), we take this LG\- 
like mode, denoted by El{t,9) for the field distribution, 
as a background potential and solve nonlinear wave equa- 
tion described in Eq. (TTJ) with an additional self-induced 
refractive index modulation Vl(v,9) = \El{t, 9)\ 2 : 



2C(I - 1) 
1 + |£| 2 



AE + d 2 ± AE + [V(r, 9) + V L (r, 9)]AE = 0.(2) 



Here we treat the electrical field E(r, 9) — El + AE and 
find the solution for AE self-consistently. The validity 
of this perturbation approach can only be applied to the 
case near the threshold condition for the reason that the 
unperturbed electrical field El (V, &) is suppressed below 
the lasing condition, but with a perturbed electrical field 
AE(r, 9) one has a chance to overcome the threshold and 
produce a lasing mode. 

The simulation results based on Eq. @ are shown 
in Fig. IHg-i), which successfully report nonlinear cres- 
cent surface modes with single-, double-, and quardruple- 
humped solutions by using the same parameters but just 
varying the injection current. As mentioned before, with 
the assistant of the small cavity in the center, which sup- 
press the supported LG^-like mode at the first. As the 
injection currents increase, this supported LG^-like mode 
is driven by a higher gain and has to turn in lasing too. 
But before this LG^-like mode to become a dominant 
lasing mode over the single-peak one, it acts as an effec- 
tive index modulation in the shape of concentric rings, 
which support stationary crescent surface modes as pre- 
dicted theoretically (23[. A one-to-one correspondence 
with experimental data in Fig. Hfd-f ) verifies our numer- 
ical investigations. 

Such a perturbed approach breaks down as the injec- 
tion current increases to a certain value. For well above 
the threshold condition, we solve the nonlinear mode 
in Eq. ([T]) directly with the supported linear mode in 
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Fig. [2c ) as an initial ansatz. The convergent mode so- 
lutions, which forms a A th -oidei cylindrically symmetric 
bound-state, are found numerically by the standard re- 
laxation method. Owing to the azimuthal instability of 
cylindrically symmetric higher-order bound states, these 
ring structures are prone to break into spots lying on the 
concentric rings, as illustrated in Fig. [3Jb). The num- 
ber of lobes is quite different to each other on distinct 
rings because the reason that modes with different az- 
imuthal indexes may have an equivalent growth rate of 
instabilities. As an example, we calculate the correspond- 
ing azimuthal instability spectrum for different azimuthal 
mode indexes in Fig.[3jc), which gives 4 possible nonlin- 
ear modes with the largest amplitude growth rates, i.e., 
I = 3,7, 11, and 14, respectively. By the superposition of 
these four nonlinear patterns illustrated in Fig. GUd-g), 
we find that not only the number of rings but also the 
number of lobes in the outer ring are in a good agree- 
ment for the case of a highly nonlinear pattern exper- 
imentally and numerically, as the comparison shown in 
Fig. [31a) and (b). In this case, one can see that the LG\- 
like mode is dominant in the radiation pattern, which 
later becomes a nonlinear mode with broken spots in the 
ring-structure as a manifold of nonlinear instability for 
soliton rings 24j . 

Before conclusion, we would like to address the com- 
parison of our introduced laser cavity to those with an 
annular Bragg resonators [29l . l30j . where the cavity modes 
are confined to propagate azimuthally within a cylindri- 



cal Bragg grating. Even though lasing modes within a 
multiple-ring structure, such as annual Bragg resonators, 
can be excited easily, but the mode competition among 
supported linear cavity modes results in the preference 
of a symmetry maintaining radiation pattern. On the 
contrary, the approach we used here seems unnatural at 
a first glance for assuming a refractive index modula- 
tion induced by the supported linear cavity modes. As 
long as one operates near the threshold condition, our 
proposed using intrinsic non-lasing cavity modes to ex- 
plain the generation of crescent surface waves works. It 
turns out that the preformed symmetry-breaking cres- 
cent waves on top of a supported linear cavity modes give 
us a satisfied comparison with experiments not only qual- 
itatively but also quantitatively. The results we showed 
theoretically and experimentally here can also be useful 
for studying a variety of microcavities with shapes from 
micro-disks, micro-rings, and micro-spheres. 

In summary, by introducing the concept to design an 
intrinsic non-lasing cavity mode as a background poten- 
tial, we fabricated a surface-structured VCSEL and re- 
ported the observation of crescent surface waves near the 
threshold lasing condition. A mode transition from sta- 
tionary single-, double-, and quardruple- humped crescent 
solitons to soliton rings are demonstrated numerically 
and experimentally. Without any counterparts in the lin- 
ear limit, the experimental observations and the simula- 
tion results provide an alternative but effective approach 
to access optical surface modes in semiconductor lasers. 
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